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RESUMEN 

Se analiza el desempeno del aislador de triple pendulo de friccion 
FPT8833/1 2-1 2/8-5 que se pensaba utilizar en la construccion del nuevo Centro de 
Investigaciones Cientificas y de Post Grado de la Universidad de Fuerzas Armadas 
ESPE, modelado como un sistema de un grado de libertad, ante la accion del sismo 
de diseno que recomienda la Norma Ecuatoriana de la Construccion, NEC-15, que 
tiene un periodo de retorno de 475 anos y ante un sismo maxi mo considerado, que 
tiene un periodo de recurrencia entre 1224 y 2190 anos, con la finalidad de 
determinar la estabilidad del sistema y los maximos desplazamientos utilizando los 
coeficientes de friccion nominales, mas bajos. 

Luego del analisis realizado se decidio incrementar la altura del aislador al 
FTP8833/1 2-1 2/8-6, que sera el que se utilice en el Proyecto. Finalmente se 
comprobo que para el nivel de cargas actuantes no es necesario colocar placas 
exteriores adicionales. 

Palabras Claves: Aislador triple pendulo. Sismos de Control. 
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ABSTRACT 

Seismic performance of the triple friction pendulum named FPT8833/12- 
12/8-5 is analyzed if used in the construction of the New Research and Graduate 
Center of Universidad de las Fuerzas Armadas ESPE. It is modeled as a system of 
one degree of freedom under the design earthquake action recommended by the 
Ecuadorian Construction Norm, NEC-15, which has a return period of 475 years, 
and the maximum considered earthquake, which has a recurrence period between 
1224 and 2190 years. The objective is to determine the stability of the system and 
the maximum displacements obtained using nominal friction coefficients. 

It was decided to increase the height of the bearing in order to reduce the 
maximum displacements. It was also found that, due to the level of acting loads, it is 
not necessary to use additional exterior plates. 



1 INTRODUCCION 

Seis de los ocho bloques estructurales del nuevo Centro de Investigaciones 
Cientfficas y de Post Grado, van a construirse con aisladores sfsmicos de triple 
pendulo de friccion FPT, por lo que se decidio hacerlo con un solo tipo de aislador 
a sabiendas de que por la variacion de las cargas en los diferentes elementos 
estructurales se podia haberlo hecho con tres tipos de aisladores, pero esto implica 
para el fabricante construir tres matrices diferentes, preparar los ensayos de las 
mismas para tres condiciones; todo esto encarece el valor de los aisladores, por lo 
que se decidio hacerlo con un solo aislador, a sabiendas de que en algunos casos 
va a estar sobre dimensionado pero en la mayor parte de los casos es el adecuado. 

Antes de indicar el aislador que se va a colocar, conviene mirar la figura 1 , 
en que se presenta un aislador FPT, con su respectiva nomenclatura. Basicamente 
el aislador esta compuesto por 4 placas de radios R t y coeficientes de roce /^; con 
alturas h h las mismas que estan indicadas en la figura 1 . 




Figura 1 Nomenclatura de un aislador FPT. Fuente: Fenz y Constantinou (2007) 
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Desempeno de aislador FPT8833/1 2-1 2/8-5 en el analisis sismico... 



El aislador FPT tiene 12 propiedades geometricas a saber: 4 de los radios 
de curvatura; 4 de las alturas y 4 de los desplazamientos. Por otra parte tiene 4 
propiedades mecanicas que son los coeficientes de friccion. 

Un caso muy comun de los aisladores es que: R ± = 7? 4 ; y que R 2 = f? 3 ; lo 
propio se tiene con las alturas: h 4 = h 4 , y h 2 = h 3 \ y con los desplazamientos: d 4 = 
d 4 , y d 2 - d 3 . Con esto las propiedades geometricas se reducen a 6. El coeficiente 
de roce de las placas exteriores tambien es igual, no asi el de las placas interiores. 
Al menos estas son las caracteristicas del aislador FPT8833/1 2-1 2/8-5 

Un modelo numerico que describe el funcionamiento del FPT, es considerar 
que este dispositivo, esta compuesto por tres pendulos, que se describen a 
continuacion. 



• Primer pendulo 

Al inicio del movimiento, se desliza la base rigida del aislador interior. Por lo 
tanto, el deslizamiento se produce en las superficies de radio R 2 y R 3 . Para las 
condiciones geometricas del aislador FPT8833/1 2-1 2/8-5, se tiene: 



L ± — R 2 + R 3 — {h 2 + ft 3 ) 



( 1 ) 




Ai= C 111 - M 2 ) Ti 



( 2 ) 

(3) 



Donde L t es la longitud del primer pendulo; T t su periodo fundamental; A t 
el desplazamiento lateral maximo; g es la aceleracion de la gravedad. 

• Segundo pendulo 

Corresponde al caso en que el aislador se desliza en las placas 1 y 3, cuyos 
radios son: R 1 ,R 3 . 



L 2 — Ri + R 3 — (hi + h 3 ) 




^2— (M4 — M 1 H 2 



(4) 

(5) 

( 6 ) 



Donde L 2 ,T 2 ,A 2 , son la longitude, periodo y desplazamiento del segundo 
pendulo. 

• Tercer pendulo 

El movimiento se da en las placas cuyos radios son: R 1 ,R 4 . Aquf se tiene: 
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L 3 — R 4 + R 4 — (// + t 4" t) 



( 7 ) 



T3 




( 8 ) 



Donde L 3 ,T 3 , es la longitude y perfodo de vibracion del tercer pendulo. 

Los tres pendulos indicados, corresponden a las tres primeras fases del 
modelo de Fenz y Constantionu (2007) en que consideran cinco fases. Las dos 
restantes que vendrian a ser, la cuarta el aislador trabaja en el segundo pendulo y 
en la quinta fase lo hace en el primer pendulo. 



• Capacidad de desplazamiento 



Para encontrar la capacidad de desplazamiento d t del aislador FPT8833/12- 
12/8-5 se debe tener en cuenta que los aisladores FPT se desplazan y rotan. 



d t — ( d 2 + d 3 ) - F ( d 1 + d 4 ) - — 



h Tl 



-‘T2 



( 9 ) 



Donde: d L son la capacidad de desplazamiento que tiene el aislador interior 
y el exterior, ver figura 1 ; L 1; l 2 son las longitudes del primer y segundo pendulo; L T1 
es la suma de los radios R 2 y R 3 ', L T2 es la suma de los radios R ± y R 4 . 

La nomenclatura utilizada por EPS (Earthquake Protection System) para 
identificar los aisladores es la siguiente: FPT R 1 0 1 /b 4 -R 2 /b 2 -H. Donde R lt R 2 son los 
radios de las placas 1 y 2; 0 1( es el diametro de las superficies concavas 1 y 4; b 4 
es el diametro del aislador interior; b 2 es el diametro del deslizador rigido; H es la 
altura interior del deslizador rigido. 

En el Anexo 1 se presenta los componentes de un aislador FTP, en la parte 
superior se observa la superficie concava de la placa 1 , que esta totalmente pulido; 
en la parte central se observa el aislador interior que tiene en su parte interna dos 
superficies concavas y una base rfgida; finalmente abajo se indica un aislador 
completamente terminado, el mismo que esta sellado para impedir que ingrese 
polvo al aislador y pueda deslizarse libremente. 



2 AISLADOR FPT8833/1 2-1 2/8-5 



Como se indico se utilizara un aislador FPT8833/1 2-1 2/8-5, que tiene la 
geometrfa indicada en la figura 2. Las dimensiones se presentan en la tabla 1; los 
coeficientes de roce en la tabla 2. EPS (201 5). 
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Desempeno de aislador FPT8833/1 2-1 2/8-5 en el analisis sismico... 




h<=4 m 



Tablal Geometriadel aislador FPT8833/1 2-1 2/8-5 



Identificacion 


Valor en 
pulgadas 


Valor en cm 


Radio de curvatura de las placas 1 y 4 
7?i = R 4 


88 


223.5 


Diametro de las placas 1 y 4, 0 ± = 0 4 


33 


83.82 


Diametro de aislador interior, b 4 


12 


30.5 


Radio de aislador interior R 2 = R 3 


12 


30.5 


Diametro de deslizador rigido b 2 


8 


20.32 


Altura de aislador interior H 


5 


12.7 


Espesor de topes exteriores t rl 
36-33 

tn = = 1-5" 


1.5 


3.81 


Espesor de topes interiores t r2 


0.5 


1.27 


Espesor de aislador interior en el centro 

tc2 


1.5 


3.81 


Espesor de aislador exterior en el centro 

tel 


2.35 


5.97 


Capacidad de desplazamiento d t = d 4 
36 — 12 — 2 * 1.5 
d x - - 10.5” 


10.5 


26.67 


Capacidad de desplazamiento d 2 = d 3 
12 - 8 - 2 * 0.5 
d 2 = = 1.5" 

2 2 


1.5 


3.81 


Altura de aislador interior h 2 = h 3 

H 5 ' 

^ 2 =^ = 0 = 2.5" 

2 2 2 


2.5 


6.35 


Altura de aislador 


12.7 


32.26 


Alturas h 4 = h 4 

hi = 2.5 + 1.5 = 4" 


4 


10.16 
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Tabla 2 Coeficientes de roce 



Propiedad 


M 2 = Ms 


Mi 


M 4 


Valor Superior UB 


0.015 


0.07 


0.08 


Valor Inferior LB 


0.008 


0.05 


0.06 


Valor Nominal 


0.01 


0.06 


0.07 



Con los valores de roce UB (Upper Bound) se determinan las fuerzas que 
actuan en el aislador, con los valores LB (Lower Bound) los desplazamientos, 
debido a que con estos valores se hallan los maximos valores de fuerzas y 
desplazamientos. 

En la tabla 3 se presenta la longitud del pendulo; el periodo y el 
desplazamiento esperado en los dos primeros pendulos. Se aprecia que el rango 
de variacion del periodo se halla entre 1.39 y 4.15 s, de tal manera que ante un 
sismo de baja magnitud se flexibiliza muy poco a la estructura pero ante sismos mas 
intensos la flexibilizacion es notable, llegando a 4.15 s. 

El desplazamiento se obtuvo en base a los coeficientes de roce nominales. 



Tabla 3 Calculo de longitud, periodo y desplazamiento en cada pendulo. 



Pendulo 


L, 


Ti 


A,- 


1 


Li = 12 + 12 — 5 
L 1 = 19" = 48.26 cm 


T 1 — 2n 

\ 

7) = 1. 


48.26 


A 1= (0.06-0.01)19 
A 1= 0.95" 


980 
39 s. 


2 


A 2 

= 88+12- (4+2.5) 
L 2 = 93.5" 

= 237.49 cm 


T 2 = 2tt 

A 

T 2 = 3. 


'237.49 


a 2 

= 0.95 

+ (0.07 - 0.06)93.5 
A 2 = 1.885" 


980 
09 s. 


3 


L 3 = 88 + 88- (4 

+ 4) 

L 3 = 168" 

= 426.72 cm 


T 3 = 2tt 

A 

T 3 = 4, 


'426.72 




980 
,15 s 



La capacidad de desplazamiento del aislador FPT8833/1 2-1 2/8-5, es: 



19 

(12 + 12 ) 



168 

(88 + 88 ) 



d t = (1.5 + 1.5) 



+ (10.5 + 10.5) 



= 22.4" = 56.89 cm 
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Desempeno de aislador FPT8833/1 2-1 2/8-5 en el analisis sismico... 



SHEAR 




Figura 3 Curva de histeresis corte desplazamiento de aislador FPT8833/1 2-1 2/8-5 

Fuente: EPS (2015) 



Un tema muy importante en el analisis sismico de estructuras es el de la 
peligrosidad sismica y mas aun en estructuras con dispositivos de control pasivo 
como son los aisladores, ya que son estos los que van a disipar la energia sismica 
que ingresa a la estructura. 



3 PELIGROSIDAD SISMICA Y ESPECTROS DE ANALISIS 



Se considera el espectro de diseno DBE, al propuesto por la Norma 
Ecuatoriana de la Construccion de 2015, NEC-15, que tiene un periodo de retorno 
de 475 anos; en la figura 4 se indica la forma general del espectro elastico. 

El proyecto se halla en Sangolquf, donde el factor de peligrosidad sismica, 
z = 0.4; se encuentra sobre un perfil de suelo clasificado como “C” (Bonifaz, 2015), 
en consecuencia los factores de sitio de aceleracion, desplazamiento y de 
comportamiento no lineal de suelo, son: F a = 1.2; F d = 1.11; F s = 1.11. Con estos 
valores se hallan los perfodos: T 0 ,T C , que definen el inicio y fin de la plataforma de 
aceleracion constante con las siguientes ecuaciones: 

Fa 
F a 



T 0 = 0.1 F s 



( 10 ) 
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T c = 0.55 F s 0 1 ) 

** a 



tCTRC 




Figura 4 Espectro elastico de la Norma Ecuatoriana de la Construccion. NEC-15 



Al reemplazar valores se halla T 0 = 0.10 s. T c = 0.56 s. Interesa encontrar la 
aceleracion espectral S aDBE en un periodo T = Is, (rango de velocidad constante) 
en la rama descendente del espectro que viene definida por: 



s a (g) = VzF a 




( 12 ) 



Las variables todavfa no definidas, son: 77 factor de region sfsmica, para la 
sierra vale: 77 = 2.48; el factor r = 1, para un perfil de suelo “C”. 



Al remplazar valores, se tiene: 

2.48 * 0.4 * 1.2 * 0.56 



->0777777 — 



T 

0.667 



(g) 



-(g) 



(13) 



En la figura 5 se indica el espectro de aceleraciones, elastico, para el sitio 
del Proyecto, hallado con el NEC-15. Para el analisis estructural, se considera que 
a partir de T = 5 s, los desplazamientos espectrales son constantes. 
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Desempeno de aislador FPT8833/1 2-1 2/8-5 en el analisis sismico... 




Figura 5 Espectro elastico para el Proyecto, para un periodo de retorno de 475 

anos. Fuente: NEC-15. 




m t rt n 

Hill 

E*cu*la Politecn cj Maaonal - InUtuto G*of r»ioo Jo** Eg red A 



TERREMOTO DEL VALLE DE LOS CHILLOS 
10 - Agosto - 1938 



Figura 6 Isosistas del sismo del 9 de agosto de 1938. Fuente: Egred (2009). 



En el ultimo siglo, el sismo de mayor intensidad que ha afectado al Valle de 
los Chillos, sitio en el que se encuentra el Proyecto, es el registrado el 9 de agosto 
de 1938, que tuvo una magnitud estimada de 7.1 y una profundidad focal menor a 
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30 km. (Barona, 2010). Muy probablemente la magnitud este sobre estimada si se 
observa las isosistas de este sismo, en la figura 6, en la escala de Mercalli 
modificada que tiene un valor maximo de XII, es de VII en la zona de mayor dano 
que fue Sangolquf y sus alrededores, cuyas coordenadas geograficas son: latitud - 
0.334 S, longitud -78.452 W. (Egred, 2009). En esta zona se encuentra el sitio del 
proyecto. 

No se ha realizado un estudio detenido de las fallas geologicas existentes 
en el Valle de los Chillos, en que se encuentra el Proyecto pero si se lo ha realizado 
de la ciudad de Quito que se encuentra aproximadamente a 15 km, y es importante 
presentar este trabajo para inferir el sismo MCE. 

Quito se halla sobre cuatro fallas geologicas inversas, denominadas: 
Puengasi, llumbisi-La Bota, Carcelen el Incay Bellavista-Catequilla, cuyas trazas y 
pianos de falla se presentan en la figura 7. La falla mas cercana al Proyecto es la 
de Puengasi que tiene un perfodo de recurrencia que se encuentra entre 1224 y 
2190 anos para que se produzca un sismo de magnitud M w = 6.4, valores 
encontrados aplicando el metodo de Gutenberg y Richter, modificado y truncado. 
Aguiar (2015), Rivas et al. (2014). 

Se encontraron espectros para el Sur de Quito, en los suelos tipo “C”, 
suponiendo que se genera un sismo M w = 6.4 en la falla de Puengasi, con una 
profundidad focal de 1 0 km., empleando tres modelos de movimientos fuertes y son 
los desarrollados por: Abrahamson et al. (2014); Campbell y Borzognia (2014) y 
Chiou y Youngs (2014). 

En la figura 8 se presentan los espectros medios hallados en perfil de suelo 
“C” en el Sur de Quito, con los modelos indicados en el parrafo anterior, para niveles 
de confianza del 50% (Ifnea continua) y 84% (Ifnea discontinua) tambien se indica 
el espectro hallado con el NEC-15. 




Fuente: Rivas et al. (2014) 

Para el perfodo de 1 segundo, los espectros hallados con un nivel de 
confianza del 84% superan la aceleracion espectral encontrada con el NEC-15 y si 
se divide la aceleracion hallada con el modelo de Campbell y Borzognia (201 4) para 
la aceleracion del NEC-1 5 se encuentra un valor de 1 .4 (0.94/0.667). 
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Desempeno de aislador FPT8833/1 2-1 2/8-5 en el analisis sfsmico... 



Considerando el periodo de recurrencia ya indicado se considera que el 
factor de 1.4 es adecuado para encontrar el sismo MCE a partir del sismo DBE, 
hasta contar con estudios similares a los realizados en la ciudad de Quito en el Valle 
de los Chillos. 




i/i 

2,60 

2.40 
2,20 
2,00 
1,80 
1,60 

1.40 
1,20 
1,00 
0,80 
0,60 
0,40 
0,20 
0,00 

0,00 



ESPECTROS DE CONTROL PARA QUITO FALLA PUENGASI 
(SECTOR SUR) SUELO TIPO C 



0-50 1,00 Perils) 2,00 2^0 3ff> 



Figura 8 Espectros encontrados en perfil de suelo tipo C, para el sector Sur de 

Quito. Fuente: Aguiar (201 5). 

Por lo tanto, la aceleracion espectral del sismo MCE, que se denomina 
S aM cE > para la zona de periodos largos, viene dada por la siguiente ecuacion. 

_ 0.94 

SaMCE ~ j (5) (14) 

Se destaca que el periodo es una definition estadistica que no tiene relation 
con el periodo de recurrencia. 



4 METODO SIMPLIFICADO LINEAL 



El Metodo simplificado lineal, para el analisis sfsmico de estructuras con 
aisladores de base FPT, propuesto por ASCE 7-1 0 considera que la superestructura 
va a desplazarse como cuerpo rfgido; que el periodo del sistema de aislacion es 
mucho mas grande que el periodo de la superestructura y que el amortiguamiento 
del sistema de aislacion es mayor que el de la estructura. Con estas consideraciones 
se puede considerar al sistema como de un grado de libertad con un modelo 
constitutivo del aislador bilineal, como el indicado en la figura 9. 

Para el diagrama bilineal indicado en la figura 9, en Aguiar et al. (2008) se 
demuestra que la rigidez efectiva k ef y el amortiguamiento equivalente /? vienen 
dados por las siguientes ecuaciones: 
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kef 



W uW 
" R + — 



(15) 



P 



9 + 



D 



(16) 



Donde W es el peso que gravita sobre el aislador; R es el radio de curvatura; 
/i es el coeficiente de friccion; D es el desplazamiento. 




Figura 9 Simplificacion del diagrama de histeresis Bilineal para el FPT. 

Ahora bien, en el metodo simplificado se tiene que ^ es el coeficiente de 
friccion en las superficies 1 o 4 (ver figura 1 ). 



R — R e/l + Ref A 


( 17 ) 


Refi ~ Ri ~ hi 


( 18 ) 



Por otra parte, el perfodo del sistema de aislacion, se halla con la siguiente 
ecuacion. 



T = 



W 

kef 9 



( 19 ) 



Los espectros se encuentran para 5% de amortiguamiento, ahora cuando 
se tienen aisladores en una estructura, el amortiguamiento p es mayor al 5%, por lo 
que se debe encontrar un factor B por el que se debe dividir el espectro del 5% para 
hallar el espectro con amortiguamiento p. En la tabla 4 se indican los valores de B, 
para diferentes valores de amortiguamiento p, recomendados por el ASCE 7-1 0. 



Tabla 4 Valores de reduccion de espectro por B. Fuente: ASCE 7-10 



P (%) 


5 


10 


20 


30 


40 


B 


1.00 


1.20 


1.50 


1.70 


1.90 
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Desempeno de aislador FPT8833/1 2-1 2/8-5 en el analisis sismico... 



Para encontrar el desplazamiento D se aplica la definicion de Pseudo 
espectro, con lo que se encuentra. 




Pero se conoce que: 



5a = 



5di 

~T~ 



(5) 



Luego: 




T 

B 



( 20 ) 



Donde 5 di es el coeficiente de la aceleracion de la gravedad para un periodo 
de 1 segundo. 

El calculo del desplazamiento se realiza en forma iteractiva. En efecto se 
inicia imponiendose un valor de desplazamiento, con el cual se obtiene la rigidez 
efectiva, el amortiguamiento equivalente y el periodo efectivo T, con este valor se 
halla el desplazamiento D en funcion de 5 di y B\ se compara el desplazamiento 
impuesto con el calculado y se continua hasta lograr una convergencia. Se 
recomienda aplicar el metodo de la biseccion para obtener mas rapido la 
convergencia. 

Lo descrito en el presente apartado, se realiza para el sismo DBE y MCE. 

5 METODO SIMPLIFICADO CON CURVA DE HISTERESIS DE 
FPT 



El analisis sismico no se lo realizo con el metodo indicado en el apartado 
anterior, se lo presento con el proposito de que se lo conozca, debido a su sencillez. 
Se lo hizo con el metodo simplificado pero trabajando con verdadera curva de 
histeresis del FPT. Fenz y Constantinou (2007, 2008), Fadi y Constantinou (2010), 
Constantinou et al. (2007), 

En realidad lo que cambia, es que en lugar de trabajar con la curva de 
histeresis bilineal indicada en el apartado anterior, se lo hace con el modelo que se 
halla en los aisladores FPT, que se describen muy rapidamente a continuacion. 

Fenz y Constantinou (2007,2008) consideran cinco fases en el 
funcionamiento del FPT, cuatro de las cuales se indican en la figura 1 0. La primera, 
segunda y tercera fase corresponden a tres pendulos indicados al comienzo de este 
articulo; en la fase cuatro el aislador interior topa uno de los topes exteriores y 
trabaja como el pendulo dos y en la fase cinco el aislador interior topa el borde 
inferior y superior, por lo que el movimiento se da en el pendulo 1 . 
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Figura 10 Fases del modelo de Fenz y Constantinou (2007) para el aislador FPT. 



Al considerar el caso frecuente de la geometrfa de los FPT en que se 
reducen las propiedades de 12 a 6, Me Vitty y Constantinou (2015) presentan un 
modelo mas sencillo de analisis en que se reduce el comportamiento de seis fases 
a tres fases. Para cada una de las fases se tienen ecuaciones para definir la ri^dl^z 
efectiva en funcion del desplazamiento esperado, que de igual manera se halla en 
forma iteractiva. Se encuentra tambien el amortiguamiento equivalente, empleando 
la siguiente ecuacion. 



Donde E D es la energia disipada en un ciclo de histeresis; E L es la energia 

elastica. 

6 DESPLAZAMIENTOS ESPERADOS 

En la tabla 1 se presentan las caracteristicas y propiedades del aislador 
FPT8833/1 2-1 2/8-5, cuyos coeficientes de friccion o de roce son los indicados en la 
tabla 2. El calculo de los desplazamientos laterales se realizo para el sismo DBE y 
MCE ya indicados, trabajando con los valores nominales del coeficiente de friccion, 
pero para que se tenga mayor desplazamiento se considero el menor valor de la 
friccion de las placas exteriores. Es decir, se calculo con los siguientes valores: 



Ai = M4 = 0-06 
Az = A s = 0-01 
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Figura 11 Desplazamiento esperado para el sismo DBE y MCE. 

Por otra parte, la carga promedio que Mega al sistema de aislacion es de 
268 T. Para esta carga los desplazamientos esperados, considerando el efecto de 
torsion que no se ha indicado en este articulo pero se lo tomo en cuenta en el calculo 
se indican en la figura 1 1 . 

Para el sismo DBE el desplazamiento lateral maximo es D - 21" = 53.3 cm, 
y para el sismo MCE es D = 22" = 55.8 cm. Son valores bastante altos, pero de eso 
se trataba ya que los coeficientes de friccion con los que se calculo son muy bajos 
para precisamente hallar desplazamientos altos pero lo importante es conocer que 
el aislador si resiste un sismo con perfodo de recurrencia entre 1224 y 2190 anos. 

Los valores indicados en la parte superior izquierda de la figura 11 
corresponden a las coordenadas del diagrama de histeresis: desplazamiento en 
pulgadas y relacion fuerza lateral para carga vertical (F/W). 

7 CALCULO DEL ESPESOR DE LA PLACA 



El objetivo del presente apartado es conocer si se necesita una placa 
exterior al aislador o si el valor de t cl = 2.35" (espesor de las placas 1 y 4 en el 
centro) es suficiente, (EPS, 201 5). En caso de no ser necesario la placa exterior se 
colocara directamente el aislador a las vigas de acero o a las columnas o pedestal 



Michel Constantinou, Roberto Aguiar, Enrique Morales, Pablo Caiza 



16 



de hormigon. El calculo se realiza en dos sitios, el primero cuando el aislador se 
encuentra en el centro (solo cargas verticales) y el segundo cuando se halla en los 
extremos (cargas verticales mas sismo). 

• Espesor de la placa en el centro 

Se aplica el metodo de la columna equivalente que consiste en considerar 
al aislador interior como una columna con dos placas de acero en sus extremos de 
ancho t eff . 



teff - 



trl + ti 



( 22 ) 



Donde t b , es el espesor de la placa exterior al aislador, que se va a 
considerar igual a cero, ya que no se necesita como se vera a continuacion; el 
significado de t cl ya fue indicado. Sobre estas dos placas gravita una carga axial P u , 
que genera una presion F bu en un circulo de diametro h t . 



Fbu 



= 0c *0.85 U 




(23) 



Donde 0 c es el factor de minoracion de esfuerzos en el hormigon que es 
igual a 0.65; f c ' es la resistencia a la compresion del hormigon; — es un factor de 

Ai 

confinamiento del hormigon que se considera igual a 2. 

El diametro h 1 se calcula con la siguiente ecuacion. 



F hl , = 



Pu _4 Pg 

nbl nbl 
~T~ 



b i = 



4 P„ 



nF, 



bu 



(24) 



Si no se halla en la forma indicada se puede asumir un valor mayor de V 
En la parte superior de la figura 1 2 se presenta el aislador en la parte central y en la 
parte inferior el area de la presion en la placa que genera un momento M u en la cara 
del aislador interior, que se obtiene con la siguiente ecuacion. 



M U = Fbu Y + Fbu ' “ 



(£-0 



r = 



bi 



(25) 



(26) 



La ecuacion para determinar el momento de flexion es aproximada ya que 
se ha considerado una carga uniforme mas una cargar triangular, como se aprecia 
a la derecha de la figura 13, en lugar de considerar la carga variable en forma 
esferica indicada a la izquierda de la figura 1 3. Por eso se debe corregir el momento 
de flexion M u por un factor de correccion que se obtiene de la figura 14, en funcion 
de la relacion b/b t . 
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Figura 12 Metodo de la columna equivalente. Fuente: Constantinou etal. (2011). 




Figura 13 Diagramas de presiones, exacto y aproximado para calcular el momento 
de flexion, en la cara del aislador interior. Fuente: Constantinou etal. (2011). 
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Figura 14 Factor de correccion del momento de flexion por la forma de calculo de 
M u . Fuente: Constantinou et al. (2011). 

Finalmente se debe verificar que la capacidad a flexion que se va a 
denominar M Y sea mayor que la demanda M a . Se denomino M Y ya que se va a 
trabajar con el momento de fluencia de las placas. 

t 2 

M y = 0 s * f y *z = 0 S * f y * (27) 

Donde f y es el Ifmite de fluencia del material; z es el momento plastico. Las 
restantes variables han sido ya indicadas. En la tabla 5 se presenta el calculo del 
espesor de la placa, en el centro. 



Tabla 5 Calculo del espesor de la placa en el punto central. 



Variable 


Nombre 


Valor 


Diametro de aislador interior 


b 


12" 


Espesor de la placa concava en el centro. 


tel 


2.35" 


Espesor de la placa exterior 


h 


0 


Espesor efectivo de la placa 

teff = Jtcl + *b 


teff 


2.35" 


Resistencia del hormigon fj = 21 MPa 


fc 


3.0 ksi 


Relacion de confinamiento en el hormigon 


A 2 /A 1 


2 


Limite de fluencia del acero f y - 276 MPa 


fy 


40 ksi 


Factor de reduccion de esfuerzos en hormigon 


0r 


0.65 
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Factor de reduccion de esfuerzos a flexion en acero 


0s 


0.90 


Presion en el hormigon 

F bu = 0 C *O.85/ C ' ^ 
\ 1 


Fbu 


3.32 ksi 


Diametro by. Considerando Py = 268 T. 

, 4 Py 

1 Jrr F bu 


by 


15.05” 


Se impone un valor mayor a by (Criterio conservador) 


by 


20.7” 


Calculo de longitud de volado r 

b 1 — b 

r ~ 2 


r 


4.4” 


Calculo de My simplificado 

r 2 /by \ r 2 

My=F buY+ F bu (-- l) y 


My 


46.5 Kip- 
in 


Relacion h/h 1 para hallar factor de correccion de Mu 


b/by 


0.58 


Factor de correccion de figura 14 


f 


0.93 


Momento ultimo exacto 


My e 


43.27 kip- 
in 


Capacidad de Momento en la fluencia 

t 2 

My = 0s* fy*Z= 0S* fy* 


My 


49.7 kip- 
in 



El momento de capacidad es mayor que el momento de demanda. Por lo 
tanto, el espesor t cl = 2.35" es suficiente y no se require placa exterior. Se destaca 
que el momento My es menor al indicado en la tabla 5, debido a que el diametro by 
es menor al valor con el cual se obtuvo el momento de demanda. 

• Espesor de la placa en el borde 

En la figura 15 se ilustra el punto en el cual se halla el aislador debido a 
cargas verticales mas un sismo muy intenso. La parte que conlleva mayor trabajo 
es determinar el diametro by, una vez encontrado este valor se procede en forma 
similar al anterior pero considerando los siguientes factores de minoracion de 
esfuerzos 0 C = 1, 0 S = 1. 

En la figura 1 6 se ilustra el proceso de calculo, a la izquierda se presenta el 
diagrama de presiones con la carga actuando en la posicion desplazada se aprecia 
que una parte del diagrama de presiones queda afuera de la placa del aislador y 
esto no puede darse, por lo que Constantinou ef a/ (2011) recomienda mantener la 
misma area de distribucion de presiones pero alargandola en el sentido contrario 
como se observa en la parte central de la figura 16, este alargamiento conduce a 
tener una distribucion de presiones en forma elfptica con eje menor a t y eje mayor 
V 
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b 




Figura 15 Posicion desplazada del aislador interior 

La determinacion del eje menor % esta dada por la geometria del aislador 
y por el desplazamiento esperado. 

-(H- (28) 

El eje mayor se encuentra con la siguiente expresion 



4 P u (29) 

&i = 

TTCLlFfou 

En lugar de calcular b x en la forma indicada se puede imponer un valor 
siempre y cuando sea mas alto que el valor que se halla con la ecuacion anotada. 
Finalmente se calcula el momento ultimo con la distribucion de presiones que se 
encuentra en el eje mas largo (area achurada de figura 16) y se procede en forma 
similar al caso anterior. 

En la tabla 6 se resume el calculo del espesor de la placa en los extremos. 
Nuevamente no se considera la existencia de placa exterior. 
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Pu 




Area de presion 
fuera de la placa 




b, 



Figura 10. Presiones en la placa ante desplazamiento horizontal. 
Fuente: Constantinou etal. (2011) 



Tabla 6 Calculo del espesor de la placa en el extremo. 



Variable 


Nombre 


Valor 


Diametro de aislador interior 


b 


12" 


Espesor de la placa concava 


tel 


2.35" 


Espesor de la placa exterior 


t b 


0 


Espesor efectivo de la placa 

t eff = Jtcl + 


teff 


2.2” 


Resistencia del hormigon f c ' = 21 MPa 


fr 


3.0 ksi 


Relacion de confinamiento en el hormigon 


A 2 /A 1 


2 


Limite de fluencia del acero f y = 276 MPa 


fy 


40 ksi 


Factor de reduccion de esfuerzos en hormigon 


0c 


1.00 


Factor de reduccion de esfuerzos a flexion en acero 


0c 


1.00 


Presion en el hormigon 

F bu =0c* 0.85 / c ' ^ 
V 1 


F bu 


5.10 ksi 


Se impone un valor de b , (Criterio conservador) 


b 1 


18” 


Calculo de longitud de volado r 

b t — b 

' ~ 2 


r 


3” 


Calculo de M u simplificado 

r 2 /fei \ r 2 

Mu = F bu -+ F bu (--l)~ 


M u 


30.6 Kip- 
in 


Relacion b/b , para hallar factor de correccion de Mu 


b/b , 


0.67 


Factor de correccion de figura 14 (Valor conservador) 


f 


1.00 


Momento ultimo exacto 


Mue 


30.6 kip- 
in 
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Capacidad de Momento en la fluencia 


My 


49.7 kip- 


t 2 

My = 0 S * fy*Z= 0s* fy* 




in 



Nuevamente la capacidad a flexion es mayor que la demanda. Por lo tanto 
no se requiere placa extra y las dimensiones del aislador son las adecuadas. 



8 COMENTARIOS Y CONCLUSIONES 

Se ha presentado las caracteristicas geometricas y mecanicas de un 
aislador FPT8833/1 2-1 2/8-5 y se ha determinado las fases en las que va a trabajar 
el aislador para un mejor entendimiento del funcionamiento de los aisladores de 
triple pendulo de friccion. 

Por otra parte se ha definido el espectro de diseno DBE para el sitio en que 
se halla el Centro de Investigaciones Cienttficas de la Universidad de Fuerzas 
Armadas ESPE, en el Valle de los Chillos en Ecuador, siguiendo lo prescrito en la 
NEC-15, que tiene un pertodo de retorno de 475 anos y se ha determinado el 
espectro para el sismo maxi mo considerado MCE, asociado a un sismo de magnitud 
M w = 6.4 con profundidad focal de 10 km, en la fall a de Puengast que es la falla 
mas cercana al Proyecto que ha sido estudiada con detenimiento. Este sismo tiene 
un pertodo de recurrencia que se encuentra entre 1224 y 2190 afios. Del estudio 
realizado se determino que el espectro para el sismo MCE es igual a 1.4 por el 
espectro del sismo DBE. 

Por otra parte se encontro el desempeno del aislador ante los espectros 
DBE y MCE, considerando los valores nominales, del coeficiente de friccion, 
entregados por EPS (2015) y se hallo que los desplazamientos laterales son 
bastante altos. Por lo que se decidio incrementar la altura del aislador para tener 
mayor rigidez y menores desplazamientos. El que se va a utilizar es FPT8833/12- 
12 / 8 - 6 . 



Final mente se comprobo que con las dimensiones del aislador FPT8833/1 2- 
12/8-5, no se requiere colocar una placa sobre los aisladores sino que deben 
acoplarse directamente a las vigas de acero que conforman la losa de aislacion. 
Para el caso en que se tiene los aisladores sobre columnas de hormigon si se 
necesita una placa exterior por el sistema constructive. 
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ANEXO 1 




PLACA INFERIOR CONCAVA DE RADIO Y DIAMETRO 0 4 




AISLADOR INTERIOR DE DIAMETRO Y COEFICIENTE DE ROCE ju 4 




AISLADOR TERMINADO 




